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Using the correlation matrix formalism wc study the temporal aspects 
of the Warsaw Stock Market evolution as represented by the WIG20 index. 
The high frequency (1 min) WIG20 recordings over the time period between 
January 2001 and October 2005 are used. The entries of the correlation ma- 
trix considered here connect different distinct periods of the stock market 
dynamics, like days or weeks. Such a methodology allows to decompose the 
price fluctuations into the orthogonal eigensignals that quantify different 
modes of the underlying dynamics. The magnitudes of the corresponding 
eigenvalues reflect the strengths of such modes. One observation made in 
this paper is that strength of the daily trend in the WIG20 dynamics sys- 
tematically decreases when going from 2001 to 2005. Another is that large 
events in the return fluctuations arc primarily associated with a few most 
collective eigensignals. 

PACS numbers: 89.20.-a, 89.65. Gh, 89.75.-k 



1. Introduction 

Nature of the temporal correlations in financial fluctuations constitutes 
one of the most fascinating issues of the contemporary physics. The pure 
Brownian-type motion [1] is definitely not an optimal reference [2] . Already 
the correlations in the financial volatility remain positive over a very long 
time horizon [3]. Even more involved are higher order correlations that give 
rise to the financial multifractality [1]. In this context a mention should 
also be given to the concept of financial log-periodicity - a phenomenon 
analogous to criticality in the discrete scale invariant version [5l [6] . 

One more approach, initiated in Ref. [7], to quantify the character of 
financial time-correlations [8] is to use a variant of the correlation matrix. 
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In this approach the entries of the corresponding matrix connect the high- 
frequency time series of returns representing different disconnected time- 
intervals hke the consecutive days or weeks. The structure of eigenvectors 
of such a matrix allows then to quantify several characteristics of time cor- 
relations that remain unvisible by more conventional methods. 
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Fig. 1. The WIG20 (Warsaw Stock Exchange), the S&P500 and the Nasdaq indices 
from 2001.01.02 until 2005.10.31. 

Using this methodology here we present a systematic study for the Pol- 
ish stock market index WIG20 over the period 02.01.2001-31.10.2005. The 
corresponding WIG20 chart, expressed both in terms of the Polish Zloty 
(PLN) and in terms of the US$, versus two world leading stock market in- 
dices: the Nasdaq and the S&:P500, is shown in Fig. 1. The Warsaw Stock 
Exchange trading time during this period was 10:00-16:10 and the WIG20 
recorded with the frequency of 1 min. 

2. Formalism 

In the present study the correlation matrix is thus defined as follows. 
To each element in a certain sequence of relatively long consecutive time- 
intervals of equal length K labeled with a one uniquely assigns a time series 
XaiU)-, where ti {i = 1, 2, K) is to be understood as discrete time counted 
from the beginning for each a. In the financial application Xa{ti) is going 
to represent the price time-series, a the consecutive trading days (or weeks) 
and ti the trading time during the day (week). As usual it is then natural 
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to define the returns Raiti) time-series as Ra{ti) = InXaiU + t) — lnxQ,(tj), 
where r is the time lag. The normahzed returns are defined by 

r^{U) = - {Ra{U))t 

V 

where v is the standard deviation of returns over the period T and = 
a'^{Ra) = {Rait))t — {Rait))t j and (. . .)t denotes averaging over time. 
One thus obtains N time series ra{ti) {a = 1,...,A^) of length T=K-1, 
i.e. an N X T matrix M. Then, the correlation matrix is defined as C = 
(l/T) MM'^. By diagonalizing C 

Cv'^ = Afcv^ (2) 

one obtains the eigenvalues (A; = 1, A^) and the corresponding eigen- 
vectors = {v^}. In the limiting case of entirely random correlations the 
density of eigenvalues pcW is known analytically [Sllinilll], and reads 

pcW = 2;^^ X ' 

Ar:=^T2(l + l/Q±2yi7Q), (3) 

with Xmin < A < Xmax, Q = T/N > 1, and where cj^ is equal to the 
variance of the time series (unity in our case) . 

For a better visualization, each eigenvector can be associated with the cor- 
responding time series of returns by the following expression: 

N 

Zk{ti) = Y,v'^ra{ti), k=l,...,N; i = l,...,T. (4) 

a=l 

These new time series thus decompose the return fluctuations into the or- 
thogonal components that reflect distinct patterns of oscillations common 
to all the time intervals labeled with a. They are therefore called the 
eigensignals [71[8]. 



3. Results 

3.1. Correlations among trading days 

The above methodology is here applied to the WIG20 1 min recordings 
during the period between January 02, 2001 and October 31, 2005. This 
whole time period is split and analysed separately for the consecutive calen- 
dar years y that it covers (3^ = {2001, 2002, 2003, 2004, 2005}). The number 
Ny of trading days correspondingly equals 249, 243, 249, 255, and 210. The 
WIG20 intraday variation is systematically taken between the trading time 
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10:01:30 (at this time the index is always already determined) and 16:10:00. 
This corresponds to T = 368 during one trading day. Using these data sets 
we construct the Ny x Ny correlation matrices C^;. 
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Fig. 2. Empirical eigenvalue spectrum of the correlations matrices Cy (vertical 
black lines) calculated for WIG20(Warsaw Stock Exchange) index over the five 
consecutive calendar years. The noise range, as determined by a random Wishart 
matrix with Q = 368 /Ny, is indicated by the shaded field. 

The structure of eigenspectrum of such matrices for all the five calen- 
dar years is shown in Fig. 2. The pure noise range - as prescribed (Eq. ^) 
by the corresponding Wishart ensemble of random matices |10l [TT] - is in- 
dicated by the shaded area. As one can see, the majority of eigenvalues of 
our empirical correlation matrices are located within this area which signals 
that noise is dominating. Typically there exist however several eigenvalues 
that stay significantly above it. They are associated with some collectivity 
effects that in the present case are to be interpreted as an appearance of 
certain repeatable structures in the intraday dynamics of financial fluctua- 
tions. Definitely one such structure is the daily trend. As far as the WIG20 
dynamics is concerned it is however even more interesting to see that when 
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going from 2001 to 2005 those large eigenvalues gradually decrease and get 
closer to the noise area. This effect is more systematically documented 
in Fig. 3 which shows the evolution of the four largest eigenvalues of the 
N xN = 250 X 250 correlation matrix, which corresponds to 250 consecutive 
trading days, and this time window is moved with a step of one month (20 
trading days). 

The structure of eigenspectrum is expected to be closely related to the 
distribution of matrix elements of the correlation matrices. For the same 
five calendar years as in Fig. 2 the corresponding distributions are shown 
in Fig. 4 versus their Gaussian best fits. Indeed, in 2001 this distribution 
deviates most from a Gaussian and even develops a power-law tail with the 
slope of 7 ~ 6 {P{x) ~ x~"'). It is in this case that the largest eigenvalue 
of the correlation matrix is repelled most (upper most left panel of Fig. 2) 
from the rest of the spectrum due to an effective reduction of the rank of the 
matrix [T^] . Later on the distribution of the matrix elements is much better 
fit by a Gaussian and some deviations remain on the level of essentially 
single entries. 
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Fig. 3. Eigenvalues evolution of the sequence of N x N = 250 x 250 WIG20 correla- 
tion matrices translated with the step of one month. The dashed line corresponds 
to noise level. 

An optimal way to visualize the character of repeatable intraday struc- 
tures is to look at the eigensignals as defined by the Eq. They can 
be calculated for all the eigenvectors. An explicit numerical verification 
confirms that they are orthogonal indeed, i.e., the correlation matrix con- 
structed out of them is diagonal. 

The most relevant examples of such eigensignals - corresponding to the 
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Fig. 4. Distribution of matrix elements C'a.a' of the Ny x Ny correlation matrices 
Cy for the WIG20 variation during the intraday trading time 10:01:30-16:10:00. 
The solid lines indicate the power law fits to the tails of C2001 with the power 
index 7 = 6 {P{x) ~ x~'^). The dashed lines corresponds to a Gaussian best fit. 



two largest eigenvalues, for all the five years considered here - are shown 
in Fig. 5. They both display a strong enhancement of market activity just 
after the opening and lasting up to 60 min. 

Interestingly, an analogous enhancement before closing as observed [8] 
for the other markets, in case of the WIG20 can be seen only rudimentary. 
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Fig. 5. Intraday eigensignals corresponding to the two largest eigenvalues (Ai, A2) 
calculated for five calender years of WIG20 (Warsaw Stock Exchange) index vari- 
ation during the intraday trading time 10:01:30-16:10:00. The last graph is the 
same intraday eigensignal but corresponding to A200 for the year 2001. 
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For comparison one typical eigensignal (2200 (ij)) corresponding to the bulk 
of eigenspectrum is shown in the bottom of Fig. 5. Its amplitude of oscil- 
lations can be seen to be about an order of magnitude smaller than for the 
previous leading eigensignals. 



3.2. Correlations among trading weeks 

Ability of the above formalism to detect and decompose some potential 
repeatable structures in the financial patterns prompts a question of correla- 
tions among different trading weeks. Our WIG20 data set (Fig. 1) comprises 
N = 207 full trading weeks and thus allows to construct a 207 x 207 cor- 
relation matrix. The lengths T of the corresponding time series of 1 min 
returns between monday opening and friday closing equals 1840. The upper 
panel of Fig. 6 shows the resulting spectrum of eigenvalues. Most of them 
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Fig. 6. Top - Empirical eigenvalue spectrum of the 207 x 207 correlation matrix 
calculated among the weekly time intervals for the whole period of the WIG20 
recordings as shown in Fig. 1. The corresponding noise range of a random Wishart 
matrix with Q — 1840/207 ( Xmax ~ 1-78 and Xmin ~ 0.44) is marked by the 
shaded field. Bottom - Intraweek (monday 10:01:30 - friday 16:10:00) eigensignal 
associated with the largest eigenvalue Ai. 
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fall into the Wishart random matrix spectrum range (shaded area) but at 
least three eigenvalues of our empirical correlation matrix stay apart. The 
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Fig. 7. Cumulative distributions of the moduli of normalised intraweek eigensignals. 
Dashed line corresponds to a Gaussian distribution and the solid line indicates a 
slope corresponding to the inverse cubic power-law. 

eigensignal associated with the largest eigenvalue (Ai) is shown in the lower 
panel of this Figure. It shows an enhanced market activity at the connec- 
tions between the days. Interestingly however this activity is much stronger 
in the middle of the week than in its beginning or in the end. 

Finally, such a decomposition of financial fluctuations allows an instruc- 
tive insight into the statistics of returns distributions. This in itself consti- 
tutes one of the central issues of econophysics. The related well identified 
stylized fact is the so-called inverse cubic power-law [131 El IS] • There exist 
also some consistency arguments that favor this law |16] . Fig. 7 shows the 
cumulative return distributions associated with several weekly eigensignals. 
Those distributions that originate from the bulk of eigenspectrum can be 
seen not to deviate much from a Gaussian even for such a short time lag 
of 1 min. The fatter tails result from the fiuctuations filtered out by the 
eigensignals connected with the largest eigenvalues. In particular, the most 
extreme events can be seen in the first eigensignal, the one whose main 
component constitute fiuctuations commonly considered a daily or a weekly 
trend. 

4. Summary 

The way of using the correlation matrix formalism, as presented here, 
opens a promissing novel window to view the character of financial fluctu- 
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ations. In particular the related concept of eigensignals allows to filter out 
all repeatable synchronous patterns in the market daily or weekly activity. 
They are connected with a few largest eigenvalues of the corresponding cor- 
relation matrix. It is those eigensignals that appear to be responsible for the 
fat tails in the return distributions. The overwhelming rest of the spectrum 
stays within the borders prescribed by the random ensemble of Wishart ma- 
trices and fluctuations of the corresponding eigensignals are essentially of 
the Gaussian type. As far as the WIG20 dynamics is concerned it is inter- 
esting to notice a gradual weakening of the daily trend effects when going 
from 2001 to 2005. A question remains whether this effect is characteristic 
to this specific market or it takes place in the other markets as well. 

REFERENCES 

[I] L. Bachelier, Ann. Sci. dc I'Ecolc Norm. Sup. Ill- 17, 21 (1900) 
[2] B. Mandelbrot, J. Business 36, 294 (1963) 

[3] R. Mantegna, H.E. Stanley, Introduction to Econophysics: Correlation and 
Complexity in Finance (Cambridge Univ. Press, Cambridge, England, 2000) 

[4] P. Oswi§cimka, J. Kwapien, S. Drozdz, Physica A 347, 626 (2005); J. Kwapien, 
P. Oswi§cimka, S. Drozdz, Physica A 350, 466 (2005) 

[5] D. Sornette, Why Stock Markets Crash (Princeton Univ. Press, Princeton, 2003) 

[6] S. Drozdz, F. Griimmer, F. Ruf, J. Speth, Physica A 324, 174 (2003); M. Bar- 
tolozzi, S. Drozdz, D.B. Leinweber, J. Speth, A.W. Thomas, Inter. J. Mod. Phys. 
C 16, 1347 (2005) 

[7] J. Kwapien, S. Drozdz, A.A. loannides, Phys. Rev. E 62 5557 (2000) 

[8] S. Drozdz, J. Kwapien, F. Griimmer, F. Ruf, J. Speth, Physica A 299, 144 
(2001); J. Kwapien, S. Drozdz, F. Griimmer, F. Ruf, J. Speth, Physica A 309, 
171 (2002) 

[9] V.A. Marchenko, L.A. Pastur, Math. USSR-Sbornik, 1, 457 (1967) 
[10] A. Edelman, SIAM J. Matrix Anal. Appl. 9 543 (1988) 

[II] A.M. Sengupta, P.P. Mitra, Phys. Rev. E 60, 3389 (1999) 
[12] S. Drozdz, M. Wojcik, Physica A 301, 291 (2001) 

[13] P. Gopikrishnan, V. Plerou, L.A.N. Amaral, M. Meyer, H.E. Stanley, Phys. 
Rev. E. 60, 5305 (1999); V. Plerou, P. Gopikrishnan, L.A.N. Amaral, M. Meyer, 
H.E. Stanley Phys. Rev. E. 60, 6519 (1999) 

[14] S. Drozdz, J. Kwapicii, F. Griimmer, F. Ruf, J. Speth, Acta Phys. Pol. B 34, 
4293 (2003) 

[15] R. Rak, S. Drozdz, J. Kwapien, Nonextensive statistical features of the Polish 

stock market fluctuations, |arXiv:phys ics/0603071, Physica A in print 
[16] X. Gabaix, P. Gopikrishnan, V. Plerou, H.E. Stanley Nature 423, 267 (2003) 



